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Abstract 

Electrons are offered a valley degree of freedom in presence of particular lattice 
structures. Manipulating valley degeneracy is the subject matter of an emerging 
field of investigation, mostly focused on charge transport in graphene [lH4] . In bulk 
bismuth, electrons are known to present a threefold valley degeneracy and a Dirac 
dispersion in each valley. Here we show that because of their huge in-plane mass 
anisotropy, a flow of Dirac electrons along the trigonal axis is extremely sensitive to 
the orientation of in-plane magnetic field. Thus, a rotatable magnetic field can be 
used as a valley valve to tune the contribution of each valley to the total conductivity. 
According to our measurements, charge conductivity by carriers of a single valley 
can exceed four-fifth of the total conductivity in a wide range of temperature and 
magnetic field. At high temperature and low magnetic field, the three valleys are 
interchangeable and the three-fold symmetry of the underlying lattice is respected. 
As the temperature lowers and/or the magnetic field increases, this symmetry is 
spontaneously lost. The latter may be an experimental manifestation of the recently 
proposed valley- nematic Fermi liquid state [9]. 
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The principal challenge in the field of "valleytronics" [T] is to lift the valley degeneracy 
of electrons in a controlled way. In the case of graphene, a number of theoretical proposi- 
tions to generate a valley-polarized fiow of electrons [H H] have emerged. They are yet to be 
experimentally realized. Valley degeneracy has also been explored in two-dimensional semi- 
conductors such as AlAs heterostructuresO |6] and the surface states of sihcon[71 [8]. These 
systems, in contrast with the case of graphene, are host to valleys presenting a parabolic 
dispersion and a significant in-plane mass anisotropy. Recently, it has been proposed that 
in presence of such an anisotropy. Coulomb interaction energetically favors a spontaneous 
imbalance in the occupation of different valleys[S]. 

In this paper, we show that in the case of semi-metallic bismuth, a valley-polarized current 
be easily generated with a well-oriented magnetic field. The Fermi surface in bismuth in- 
cludes three cigar-shaped electron valleys lying almost perpendicular to the high-symmetry 
axis known as the trigonal axis [10]. The in- plane mass anisotropy of each valley exceeds 
200 [TT]. This exceptional feature is a consequence of a Dirac dispersion [12], with a quadratic 
correction [13] pulling drastically down the effective mass along two out of the three orien- 
tations. At the bottom of the band, experiment yields mi = O.OOllme, m2 = 0.26me and 
ms = 0.0044me[l0] in fair agreement with the tight-biding model|llj. In the last few years, 
the electronic properties of bismuth in presence of a strong magnetic has become a focus of 
attention [T5l - [T9] . The high-field phase diagram of the Dirac electrons in the extreme quan- 
tum limit [T7] presents more lines than what is the expected in the non-interacting theoretical 
picture [IHl EH]. The origin of these additional field scales is yet to be understood. 

We report on a study of angular-dependent magnetoresistance on bismuth single crystals 
with a mobility in the range of 10^ — 10''^ cm'^V^^s^^ (see the SI section for details on sample 
mobility), performed in a configuration designed to tune the contribution of each of the 
three valleys to charge transport. In a wide range of temperature and magnetic field, the 
magnetoconductivity of electrons can be quantitatively described as the sum of three con- 
tributions interchangeable by a 120 degrees rotation. As the field rotates, the contribution 
of each valley to the total conductivity can be tuned and a valley-polarized conductivity 
exceeding 80 percent can be achieved. At high temperatures and low magnetic fields, this 
simple description based on the assumption that there is no difference between the zero-field 
population of the three valleys is successful. However, as the temperature is lowered or the 
magnetic field is increased, magnetoresistance does not display the threefold symmetry of 



the underlying lattice and the conductivity cannot be described as the sum of three equiva- 
lent, rotationally-symmetrical channels. This spontaneous loss of valley degeneracy may be 
a consequence of electron-electron interaction in presence of anisotropic electron dispersion 
as recently suggested|9j. Fig. la shows the experimental configuration. Resistivity was 
measured along the trigonal axis, with a magnetic field rotating in the (binary, bisectrix) 
plane. In this configuration, the vector product of the two vectors, magnetic field, B, and 
current density, J, remains constant during the rotation. Charge carriers traveling along the 
trigonal axis are either holes or valley degenerate electrons. The response of these carriers 
to the applied magnetic field is determined by their mobility in the plane perpendicular to 
the orientation of the magnetic field. Because of the extreme anisotropy of the mobility 
tensor of each of the three Dirac valleys (Fig. lb), angular oscillations in magnetoresistance 
emerge at temperatures as high as room temperature and magnetic fields as low as 0.7 T 
(Fig. Ic). 

Fig. Id presents a typical set of data for one of the samples at B= 0.5 T at different 
temperatures. Because of the lightness of electrons uJcT is larger than unity in the entire 
temperature range. Here, ojc = ^ is the cyclotron frequency, r the scattering time, e the 
electron charge and is the effective mass in the plane perpendicular to the magnetic 
field. The dependence of the magnetoresistance on both temperature and the orientation 
of the magnetic field can be qualitatively understood by considering the magnitude of UcT. 
As seen in the figure, resistivity increases with decreasing temperature. This insulating- 
like behavior is visible for all orientations of the magnetic field[22]. This is because the 
temperature dependence of magnetoresistance is not set by the thermal evolution of the zero- 
field resistivity, which is metallic [21] (See the SI section for details) but by the magnitude of 
UcT. Since the scattering time increases with decreasing temperature, the magnetoresistance 
follows. The n/3 periodicity of the angular oscillations is also set by UcT. Magnetoresistance 
peaks each time the magnetic field is oriented parallel or anti-parallel to a bisectrix axis. 
In this configuration, there is a set of carriers traveling in the plane perpendicular to the 
magnetic field with minimal cyclotron mass (and maximal Uc)- These carriers are those 
which give rise to the largest large magnetoresistance. As the in-plane mass of the hole- 
like carriers is isotropic, one does not expect them to be affected by the orientation of the 
in-plane magnetic field. 

Down to 20 K, these angular oscillations present an equal amplitude within our exper- 
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FIG. 1: Fig. 1 The experimental configuration, the structure of the Fermi surface and 
angle-dependent magnetoresistance. a: Magnetoresistance of bismuth crystals was measured 
along the trigonal axis while the magnetic field, B, was rotating in the (binary, bisectrix) plane. 
The two vectors B and J were kept normal to each other during the rotation, b: The orientation 
of the three electron valleys with respect to the crystal axes in bismuth, c: Angular dependent 
magnetoresistance at room temperature at different magnetic fields. Oscillations are detectable 
when the magnetic field exceeds 0.7 T d: The evolution of angle-dependent magnetoresistance at 
B=0.5 T with temperature. Resistivity oscillates with vr/S periodicity. It peaks each time the field 
is parallel or anti-parallel to one of the three bisectrix axes. 



imental uncertainty. Below this temperature, an inequality in the amplitude of angular 
oscillations grows steadily as the temperature decreases. We first focus on the regime where 
the threefold symmetry is preserved. Fig. 2a presents a polar plot of the resistivity data 
at a fixed temperature and magnetic field (T= 40 K and B= 0.5 T). Such a presentation 
was first used by Mase and co-workers [23] and is particularly instructive to compare the 
rotational symmetries of the magnetoresistance and the underlying lattice. 

In a multi-valley system in presence of a magnetic field, the electric conductivity, a, is 
expected to be the sum of the contributions by individual valleys pi]. In our particular case, 
one can write: 



where the conductivity of electron pockets and the hole pocket are indexed by ei^s and 
h. The conductivity tensor, cf, is the inverse of the resistivity tensor, p and, in general, the 
precise determination of (Jzz implies knowledge of all off-diagonal resistivity components. 
However, in the particular case of bismuth, a compensated semi-metal, the off-diagonal 
components are known to be small (i.e. p^z ^ pzz and pzy ^ Pzz [21] )• Therefore, setting 
azz = P^z is approximation which is valid at least up to a few percent. 

Fig. 2b presents a polar presentation of the angular dependence of conductivity extracted 
from the resistivity. We found that the simplest fit to the experimental data is given by the 
following function with three adjustable parameters: 



We assume an angular- independent conductivity for the hole pocket a . The contribution 
of each of the three electron pockets is described by the same angular function rotated by 
27r/3. The magnetoresistance of an electron ellipsoid is lowest when the field is oriented 
along the binary axis which lies perpendicular to the longer axis of that particular ellipsoid. 
In this configuration, the magnetoconducivity of the ellipsoid in question acquires its largest 
value, namely cxfem. As the field rotates, magnetoresistance increases and the conductivity 
decreases until the field becomes parallel to the longer axis of the ellipsoid, which corresponds 
to a bisectrix orientation. In this configuration, the magnetoresistance attains its peak and 
the magnetoconductance its minimal value {abis = f^)- Thus, the parameter r is a measure 
of the anisotropy of magnetoconductivity of an ellipsoid at a given temperature and field. 




(1) 



i=l-3 





In a first approximation, the magnetoresistance is set by the cyclotron mass of electrons 
for a particular orientation of the magnetic field. Now, the cyclotron mass along the binary 
axis (m^*" = ^Jrfv^n^^ is 15 times larger than along the bisectrix axis(m^** = y/m^m^). 
Thus, Eq. 2 assumes an anisotropy which matches the known structure of an electron 
ellipsoid. 

Fig. 2b displays in a polar representation the best fit to the data by this equation. As 
seen in the figure, the fit is satisfactory and allows a quantified observation about the relative 
contribution of each pocket to the total conductivity. Fig. 3a shows how the magnetic field 
tunes the contribution of different valleys at T=40K and B=0.5 T. Since the contribution 
of the hole pocket is a small fraction of the total conductivity, for all orientations, charge 
conductivity is dominated by the flow of electrons and not holes. As the field rotates, the 
contribution by different valleys are modulated. When it is oriented along a binary axis, one 
out of the three electron pockets dominate the total conductivity. As seen in the figure, the 
valley-polarized conductivity exceeds 80 percent in this configuration. This is a consequence 
of the large magnitude of r, in turn a consequence of the large in-plane anisotropy of the 
cyclotron mass for each valley in bismuth. In a wide range of temperature and magnetic 
field, when the threefold symmetry of the lattice is preserved, Eq. 2 was found to be a 
successful fit to the data(See the supplementary Information section). This allowed us to 
extract the three parameters a^, abin and r for each temperature and magnetic field. Fig. 
3b plots the temperature dependence of r for different magnetic fields. As seen in the figure, 
the peak value for r (i.e. ~ 40) is attained in the intermediate temperature range centered 
around 40 K. This is where the in-plane magnetic field is most efficient for generating a 
valley-polarized current. 

Below this intermediate temperature range, r decreases as the temperature lowers. This 
suggests a reduction in the mobility anisotropy below 30 K. Mobility is set by the ratio 
of the scattering time to the effective mass. The scattering time can be anisotropic in its 
own turn and this anisotropy can be temperature dependent. In a first approximation, its 
anisotropy damps the anisotropy of the effective mass. When the field is along bisectrix, the 
effective mass in the perpendicular plane is lighter but since the phase space for scattering 
is reduced, the scattering time is longer. Below 20 K, the enhancement in the anisotropy 
of the scattering time pulls down the anisotropy of mobility and r. This would indeed be 
the case, if electrons residing in each valley are scattered by those residing in another valley 
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FIG. 2: Fig. 2 The contribution of different components of the Fermi surface to the 
total conductivity in polar coordinates, a: Polar representation of the angular dependence 
of resistivity at B=0.5 T and T= 40 K. b: Polar plot of charge conductivity (red line) according 
to the experimental data compared to a three-parameter fit (black line) following Eq. 2. The 
contribution of each electron pocket (el-e3) as well as the hole pocket (h) to the total conductivity 
is also plotted. 



through Coulomb interaction as initially suggested by Hartman|21j. 

At high temperatures, r also decreases. The oscillations gradually smear out as the 
thermal energy becomes comparable to the Fermi energy ( ~ 22meV) for the electrons pUl 
[TT] . Remarkably however, in this temperature regime, a large r is gradually restored as 
the magnetic field increases. This suggests that orbital magnetoresistance can become an 
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FIG. 3: Fig. 3 Field-induced valley-polarization across a wide temperature range, a: 

Angular dependence of the contribution of each electron valley as well as their sum, Yl'^e' 
total conductivity, Ototai- When the field is along a binary direction, because of the large r, one of 
the valleys easily dominates the total conductivity. Moreover, the total conductivity by electrons is 
much larger than the contribution of the holes, b: The temperature dependence of r for different 
magnetic fields. A large r is attained in the intermediate temperature range. For each temperature, 
r is obtained by finding the best fit between the data and Eq. 2. Error bars represent a window 
defined by the standard deviation between the fit and the data. 

effective valley valve, even when the system exits the degenerate Fermi regime, provided 
that the magnetic length(£B = \J~^) becomes shorter than the thermal de Broglie length 
(A = ^2iTm* ksf '^ ' '^^^^ interesting regime has been poorly studied and deserves further 
exploration. 

The findings reported above show that, in presence of anisotropic dispersion, the orbital 
magnetoresistance can tune the contribution of different valleys to the total conductivity. 




In multi- valley systems with modest mass anisotropy such as Si (111) or AlAs or on the 
(111) surface of bismuth itself |25] . this effect can be used to inject a valley-polarized flow to 
the two-dimensional system. However, because of the small mass anisotropy (^ is ~ 5 in 
AlAs[.5j and 3 in Si(lll)|,7J), the effect is expected to be much less drastic than in bulk 
bismuth. On the other hand, due to the absence of in-plane anisotropy no such effect is 
expected in graphite or multi-layer graphene. 

We now turn to the low-temperature high-fleld regime. Fig. 4b shows the thermal evo- 
lution of normalized Pzz{4')- As the temperature is lowered, the c3 symmetry, the threefold 
symmetry of the underlying lattice, is clearly lost. We have repeated the same kind of mea- 
surements in six different crystals and with two different set-ups. As seen in the SI section, 
each time, the threefold symmetry, clearly present in the high-temperature data was absent 
in the low-temperature data, however the pattern emerging at low temperature differed from 
one sample to the other. 

A similar evolution can be seen with increasing magnetic fleld at a flxed temperature as 
seen in Fig. 4a, which shows the data obtained on another sample. The temperature was 
kept at 10 K and the magnitude of the rotating magnetic fleld was steadily increased. As 
seen in the flgure, the c3 symmetry is preserved at B= 0.02 T and lost at B= 2 T. 

We also carried out a two-axis rotation experiment on yet another sample. In this ex- 
periment, one could rotate the sample in two perpendicular planes. The relative orientation 
of the crystal axes and the magnetic fleld was scanned both in the (binary, bisectrix) plane 
and in out of it. The results are presented in Fig. 5. As seen in the flgure, the threefold 
symmetry clearly present at 30 K is lost when the single crystal is cooled down to 1.5 K. 
This set of data clearly rules out the possibility that the loss of threefold symmetry is due 
to a residual misalignment. 

Because of the lightness of electrons, the cyclotron energy in bismuth becomes large 
even with a modest magnetic fleld. The Shubnikov- de Haas effect (quantum oscillations 
of resistivity, which are periodic in B^^) are visible in our data. But our flndings here 
primarily concern the large non-oscillating background (See the supplementary section for 
details). In other words, in the panels of Fig. 4, Landau quantiflcation leads to the flne 
structure emerging near each lobe at low temperature and/or high fleld. However, the lobes 
themselves are not caused by the passage from one Landau level to another. 

A possible explanation for the loss of threefold symmetry may be the nematic valley 




FIG. 4: Fig. 4 Spontaneous loss of threefold symmetry by lowering the temperature 
or increasing the magnetic field, a: Polar representation of the angular dependence of the 
magnetoresistance at a fixed temperature (T= 10 K) for different fields. At low fields the three- 
fold symmetry of the lattice combined to the inversion symmetry generates a sixfold rotational 
symmetry. As the magnetic field increases, Pzzi<P) loses the three- fold rotational symmetry of 
the underlying lattice, b: Same at a fixed field (B= 0.5 T) for different temperatures. As the 
temperature decreases the three- fold rotational symmetry is lost. 

scenario according to which Coulomb interaction favors an unequal occupation of valleys in 
presence of anisotropic dispersion [9]. A nematic liquid of electrons is formed when a metal 
spontaneously looses its rotational symmetry while keeping its translational invariance|26]. 
Nematic Fermi liquids have been reported in strongly correlated electron systems such as 
quantum Hall systems [27]. Sr'^Ru907P5] and more recently in URuoSioPS]. 

In our data, we do not detect any sharp anomaly in the temperature dependence at a given 
magnetic field and therefore, it is very hard to extract a critical temperature. Increasing the 
magnetic field and/or lowering the temperature regime leads to a loss of threefold symmetry 
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FIG. 5: Fig. 5. Magnetoresistance data with rotation in two perpendicular planes. 

a: Color plot of ^'^^''^^ at T=30 K and B= 0.5 T. Here 9 is the angle between the field and the 

pmax 

trigonal axis and (p is the angle between its projection and a bisectrix axis, b: Stereographic 
representation of the same data. c: Color plot of the data for T= 1.5 K and B =0.5 T.d: 
Stereographic representation of the low-temperature data. The three valleys cease to be equivalent 
at low temperatures. 

in the magnetoresistance data, which is well above our experimental margin. As seen in Fig. 
6, this loss of threefold symmetry becomes detectable below a characteristic temperature, 
which shifts upward with increasing magnetic field in agreement with what is theoretically 
expected. We did not observe any hysteresis and each crystal reproduced the same low- 
temperature pattern when the measurements were repeated (See supplementary information 
for details). This suggests that even if Coulomb interaction plays a role, the orientation 
along which the symmetry breaks is set by the crystalline imperfections of a given sample. 
Charge transport in bismuth is dominated by electron-electron scattering [21] residing in 
extremely eccentric valleys and these features appear to be a fertile ground for the emergence 
of electronic liquid crystals[2B]. 
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FIG. 6: Fig. 6. The emerging magnetoresistance anisotropy at different magnetic 
fields. 1^ as a function of temperature. Here pi refers to magnetoresistance resistivity for a field 
applied along each of the three bisectrix axes. The average resistivity is pa?) = ^ / ^ipi + P2 + Ps) and 
the valley-indexed excess of resistivity is defined as 5pi = pi — pav At high-temperature, because 
of residual misalignment between magnetic field and the trigonal plane, there is a temperature- 
independent difference between the three pi. Below, a field-dependent threshold temperature, 
this difference steadily grows and saturates to a finite value due to an intrinsic loss of threefold 
symmetry 
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Supplementary material for: 

Field-induced polarisation of Dirac valleys in bismuth 



Methods : Angular-dependent resistivity was measured both with a PPMS and with a home-made 
set-up using a piezoelectric rotator to rotate the sample and a pair of Hall sensors to determine the 
angle. Both yielded similar results. Bismuth single crystals used in this study had a typical size of 1-2 
X 1-2 X 4-5 mm^ and a "Residual Resistivity Ratio" (RRR= R(300K)/R(4.2 K)) in the range of 150 to 170. 
They were characterized as detailed below. 

1. The temperature dependence of the zero-field resistivity 




T ' ' (K' ') 

Fig. SI- Top: Out-of-plane resistivity (i.e. for a current applied the trigonal axis) as a function of in a bismuth single 
crystal used in this study. The slope of the red solid line, a guide for eyes, is 12 nficm K"^. Bottom: Same data plotted as a 
function of T^^. 

The resistivity of the sample used to obtain the data of figures 1 and 2 is 0.9 |aQcm at 4.2 K. Give the 
extremely low carrier density in bismuth (ne=nh = 3 lO" cm"^ [SI]), such a conductivity implies an 
average carrier mobility as large as <|a> =1.1 10^ cmV^ s"\ This average mobility corresponds to an 
average mean-free-path of about 0.1 mm. 

The temperature dependence of the resistivity is shown in Fig. SI and is close to T^. The slope of this 
quadratic temperature dependence (12 nQcm K'^ ) is close to what was found in previous studies[S2- 
S3]. A deviation from the T^ dependence is, however, detectable at low temperature. As seen in the 
bottom panel of the figure an exponent somewhat larger than 2 yields a better fit to the data. 
The origin of inelastic scattering in bismuth has been a subject of controversy and both electron- 
electron scattering [S2] and phonon-electron [S4] scattering have been put forward. Conventional 



phonon scattering is expected to yield a temperature dependence in temperatures below 0d/5 
(0D=12O K is the Debye temperature [S5]). This is not the case of bismuth, where the exponent of the 
power law remains close to 2 down to 0.1 K [S3]. The deviation from T^ dependence in bismuth is yet 
to be understood given the fact that in semi-metallic graphite resistivity unambiguously displays a T^ 
behavior[S6]. 



2. Low field-magnetoresistance 

In bismuth, the weak-field regime (cOc t <1) characterized by a magnetoresistance is restricted to 
very low fields[S3,S7]. Fig. S2 shows the low-field resistivity data on one of our samples. At T= 4.2 K, 
the behavior is detectable at very low fields (B < 0.4 mT). 
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Fig. S2- Low-field resistivity of another bismuth single crystal along the trigonal axis as a function of the square of the 
magnetic field. 

The magnitude of the B^ term and the range of validity both scale with mobility. As expected, the 
mobility lowers at higher temperatures and as a consequence the slope decreases and the range of 
validity of the quadratic field dependence expands. At 4.2 K, the magnitude of B^ term (0.022 Qm T^) 
yields an average mobility of 0.7 10'' cm^V"' s'^ comparable to the one estimated from the absolute 
magnitude of the resistivity. 
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FIG. S3 Top : Magnetoresistance when the magnetic field is oriented along two nominally equivalent binary axes. Bottom : 
Magnetoresistance when the magnetic field is oriented along a bisectrix and a binary axis. Quantum oscillations are visible 
on top of a monotonous background. They represent a modest change overall resistive response. The orientation of 
magnetic field has a drastic effect on the magnitude of the background. 



Angular magnetoresistance and quantum oscillations 

At low temperatures, we can resolve the quantum oscillations of the resistivity (the Shubnikov-de 
Haas effect). These oscillations are intimately linked to the topology of the pockets of the Fermi 
surface and their periodicity is affected by the orientation of the magnetic field. However, the angle- 
dependent magnetoresistance we study here is about the dependence of the non-oscillating back 
ground on the orientation of the magnetic field. This can be seen in Fig. S3. Both the bisectrix-binary 
anisotropy and the binary-binary anisotropy (i.e. the loss of the threefold symmetry at high fields 
and/or low temperatures) are associated with a drastic variation of the background resistive 
response and not with the amplitude of the quantum oscillations. 
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FIG. S4 Top: Angular-dependent resistivity in six different single crystals of bismuth at low temperatures. In all samples, the 
threefold symmetry of the lattice is lost at low temperatures. Bottom: Same at 40K. In all samples, the threefold symmetry 
is recovered, within experimental margin. 



4. Loss of threefold symmetry in different samples 

We studied six different bismuth samples cut from two larger single crystals. As seen in the top panel 
of Fig. S4, all crystals showed a departure from threefold symmetry, significantly larger than our 
experimental margin of error. The latter can be seen in the bottom panel, which shows the data for 
the same samples at 40 K. Samples a, b and c were cut from one large single crystal and samples e, f 
and g were cut from another large single crystal, with comparable residual resistivities. As seen in 
the figure, the low-temperature patterns of anisotropy differ in the two single crystals, but are 
similar in samples cut from the same single crystals. We found that the anisotropic patterns are 
robust and do not display any significant evolution with thermal recycling. Samples d and e were 
used for extensive measurements of field dependence. 

5. Loss of threefold symmetry in field scans 

Fig. S5 presents the ratio of magnetoresistance along two interchangeable bisectrix axes as a 
function of magnetic field at different temperatures in one of the crystals. As seen in the figure, at 
low temperatures, a difference appears and gradually grows with increasing magnetic field. No sharp 
anomaly is visible. However, as the temperature is decreased, the transition becomes sharper and 
one can clearly define a field above which the anisotropy (which was found to be between 7 and 15 
percent in different samples) ceases to increase further. 
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FIG. S5 The evolution of the ratio of resisitivities for field oriented along two distinct bisectrix axes with temperature. At low 
temperature the ratio saturates above a threshold magnetic field. 



6. The origin of empirical fit and its limit of applicability 

We analyzed the experimental data by a fit to equation 2. According to this equation, the 
contribution of each electron valley to the total conductivity has the same angular dependence 
rotated by 27r/3 radian. For each valley, the angular dependence is: 

a\z = abin/(l+rcos^(|)) 

This empirical fit is inspired by the well-known semi-classical expression for magnetoconductivity: 

o(B)-ao/[l+(cOcT)'] 



The angular dependence of magnetoconductivity is set by the angular dependence of the cyclotron 
frequency (governed by the anisotropy of the cyclotron mass) and the scattering time. The cyclotron 
mass is largest when the field is parallel to the binary axis and lowest when the field is parallel to the 
bisectrix axis. As a consequence, the magnetoconductivity is largest (i.e. the magnetoresistivity is 
lowest) when the field is parallel to the binary axis and is lowest when the field is parallel to the 
bisectrix axis. As seen in Fig. S6, this empirical fit is quite successful at high temperatures, but fails 
when the threefold symmetry is lost. 




FIG. S6 A comparison of experimental data of angle-dependent conductivity (red line) and the equation used to fit the data 
(black line). As seen in the data, the fit, which is quite successful at high temperature, fails below 10 K. 
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